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Abstract 

We describe a new method which allows one to evaluate the false vacuum decay rate for 
a general potential which may depend on an arbitrary number of scalar fields. 
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1 Introduction 



The problem of calculating the transition probability between two inequivalent vacua in 
quantum field theory arises in different areas of high-energy physics, cosmology, and con- 
densed matter physics. The semiclassical calculation of the false vacuum decay width was 
done in Refs. |], @, f| for the case of a single scalar field, (p(x). The corresponding path 
integral is dominated by the field configuration 4>(x) called the "bounce" and can be eval- 
uated using the saddle point method j2L [|]. The bounce, being the stationary point of the 
Euclidean action, is the non-trivial solution of the corresponding Euler-Lagrange equation 
which obeys certain boundary conditions. 

Then the transition probability per unit volume, in the semiclassical limit is 



r/v = hsy>}fe- s[ ^ n 

h 2 



det>[-di + u»m 



-1/2 



x (i + 0(h)) 



det[-d£ + U"(0)\ 

where S[<f>] is the Euclidean action of the bounce, and det' stands for the determinant with 
the zero eigenvalues omitted. 

However, a number of practical applications require a generalization of this method to 
the case that involves several scalar fields. Although, from a theoretical point of view, such a 
generalization is rather straightforward^, in practice it is impossible to find the bounce except 
in trivial cases. The problem is that the bounce is an unstable solution of some system of 
non-linear differential equations which cannot be solved analytically. Small changes in the 
initial conditions lead to significant changes in the shape of the solution. In the case of 
a single scalar field, <p(x), one can find the bounce numerically by imposing the boundary 
condition on the value of the so called "escape point", e = 0(0), and constructing the 
"overshoot" and the "undershoot" solutions which envelope the bounce. Due to the special 
topology of the one-dimensional problem, one can then constrain the bounce to any given 
accuracy by varying the "escape" point, the value of the field at the center of the bounce, 
so that the "overshoot" and the "undershoot" solutions converge. If, however, the scalar 

1 If the potential has a continuum of degenerate minima, then the determinant in equation (Q) acquires 
some additional zero modes, so that the right-hand side of (|l|) is, in fact, divergent. The false vaccum decay 
rate in this case was calculated in Ref. [Hi 
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field has n > 1 components, this approach fails because an arbitrary path connecting the (j) e 
of the "overshoot" with that of the "undershoot" does not have to go through the "true" 
escape point. 

The bounce 4>(x) is defined as a stationary point of the action S[<f>] and is the solution 
of the variational equation 5 S[<f>] = 0. However, it is not a minimum of S[<p} but rather a 
saddle point Q. 

If one could find some functional S[<f>] for which the bounce would be a minimum, then 
one could define the the function <p{x) on a lattice and minimize S[<f>] with respect to random 
variations of <fi(x). This is the basic idea of our method. We will show that by adding some 
auxiliary terms to the action one can turn the saddle point into a true minimum. First, we 
derive some useful identities which help identify the terms that may be added to the action. 
And then we prove that the saddle point of the original action is a minimum of the improved 
action. 

2 General properties of the bounce 

Let us consider a quantum field theory with the scalar potential potential U(4>i, <p n ) which 
has a local minimum at 0i = 02 = ••• = 4>n — 0, U(0,...,0) = as well as at least one 
additional (local, or global) minimum at 0j = </>?, i — 1,2, ...,n; U((j)\ , 0£) < 0. In 
the semiclassical limit |], |2|, [|, the transition probability is proportional to exp(—S[4>(x)]), 
where <j)(x) = (0i(x), is the so called "bounce", a nontrivial 0(4)-symmetric @ 
field configuration determined by the system of equations 

with the following boundary conditions: 

2 The method proposed in Ref. |g| relies on the minimization of the action on the lattice with respect 
to the variations of the bounce profile, while maximizing the same action with respect to one particular 
variation, namely the scale transformation. Since the scaling corresponds to the direction of the instability 
of the bounce []J , one can argue that it is desirable to minimize the action with respect to the variations 
orthogonal to scaling, and at the same time maximize it with respect to scaling. However, in practice, it 
is nearly impossible to identify the variations orthogonal to that related to scaling. And, in general, an 
iterative procedure of this kind does not converge to any limiting point. 
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' (d/dr)&(r)\ r=0 = 

< _ (3) 
t (fii(oo) = 

As was explained in the introduction, the chaotic instability of the solution of (0-0) with 
respect to small perturbations makes it impossible to solve the equation numerically. We 
therefore turn to the variational equation (for which (|2|) is the consecutive Euler-Langrange 
equation) : 

5S = (4) 
where S is the Euclidean action for a scalar field 4>(x): 



S = T + V (5) 

V = J U((j)(x))d 4 x (7) 

For any renormalizable theory^, the potential term in the Euclidean action V = J U(x)d 4 x 
can be written as 



V = V 2 + V 3 + V 4 (8) 

V n [<j>\ = J U n (<p(x))d 4 x (9) 

where U n {(f) is n-linear in 0. In other words, 

V n [\<f>] = J U n {\<P{x))d A x = J A n f/ n (0(x))rf 4 x = X n V n [(f)] (10) 

3 In the presence of some higher-order (non-renormilizable) terms, equations ( p3|) and (^) are modified 
accordingly to include the monomials of higher powers. However, the same method can be used to construct 
the improved action. 
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Here we used the fact that = is the local minimum and therefore no term, linear in 0, is 
present in the potential. 

By virtue of (|]), the solution <f)(x) is the stationary point of the action. Therefore, for 
any family of functions 4>\(x) such that (j)\{x)\ x=1 = (f>(x), the corresponding action S[4>\(x)} 
should have a zero derivative with respect to A at A = 1: 

^ S IM*)} L =1 = o (ii) 

In particular, one can choose 

0( A )( X ) = A p 0(A 9 x), (12) 

where p and q are arbitrary real numbers. 
Then 

4 

S[<j) w (x)] = A 2(p ~ 9) T + A pn " 49 K (13) 

n=2 

The equation (|TTj) implies 

4 

2[p-q]T+Y / [np-4:q}V n = 0^ 

n=2 

p{2T + 2V 2 + 3V 3 + W 4 ) - 2q{T + 2V) = (14) 
The relation ( |i~4|) can be satisfied for all p and g if and only if 



T + 2V = (15) 
2^ + ^ = (16) 

The first of these two identities fllS]) is well known and can be derived |], 0] using the variation 
of the type ([12]) with p = 0, g = 1. 
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Another useful identity can be obtained by considering the following variation about the 
bounce: 



(j)j(x) ■ 

4>k(x) 

It is convenient to write ([L7]) as 



lA!, .r) 
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4>i(x) - \4>j{x) 

4>j(x) + \<j)i(x) 

-(j)k(x), k=£i,j 



(17) 



z(x) = 4>i{x) + i<f>j(x) — > 2 (A) (x) = e iA (0j(x) + i(f>j(x)) 
Then it is clear that the kinetic term in the Euclidean action 



T 



u.l,- dx a 



Jd*x £ 



<d(j)k{x a ) 2 d s , d 



a dx a dx a 



(19) 



is invariant under the transformation fll8|), and therefore under the transformation fll7l) . 

It is also clear that the transformation (O) preserves the boundary conditions (Q). In 
other words, if the constraints (|3|) are satisfied for A = 0, then they are also satisfied for any 
non-zero value of A. 

We now apply the transformation fll7|) to the Euclidean action and require that it have 
zero derivative with respect to the parameter A at A = 0: 



= (d/d\)S{<p x } x=0 = {d/d\){T + V) x=0 = (d/d\)V\ x=0 



d d 
d x(4>j{x)— - 0*0)— )U(<l>)\^ x)= ^ x) 

Thus for an arbitrary pair of fa and 0,- we obtain the constraint: 



(20) 





d d 


/ d 4 x det 


drj>i dcpj 




<j>i 4>j 



U(<f>i, ...,(pr. 



(21) 



There are (n — 1) linearly-independent constraints of the form fl2~l|) for {n — 1) pairs of 
components of the n-component field <fi. 



A 2 


= 2V 2 [0(x)]4 


-V 3 [4 


>(*)] 
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= y c? 4 x det 


9<t>i 
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3 Improved Action method 

In the preceeding section we proved that the following quantities vanish identically: 

Aj = r[0(x)]+2V[0(x)] (22) 

(23) 
(24) 

In order to determine the bounce numerically, we would like to find some functional S[4>] 
for which the bounce <f>(x) would be a minimum. Then one can discretize the problem by 
defining the field <f)(x) on a lattice. The improved action, S[<p] can then be minimized with 
respect to random perturbations of the field <fi until the limiting bounce is recovered to a 
required accuracy^. 

We define S[(p} as the action S[<f>] plus some auxiliary terms which would vanish as 4>(x) 
approaches 4>(x): 

S[<f>] := S[<f>] + SKnlAi^ + « 2 , n |A 2 |P" + a 3 ^n|Af |^), (25) 

n 

where p n , n = 1,2, ... are some positive numbers, and ot n > are some arbitrary Lagrange 
multipliers. 

In particular, we will show that it suffices to take 

S[<f>] ^SM + ailAil (26) 



4 Clearly, it is not even necessary that S[<f>] be differentiable with respect to at = <f>{x). 
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but in practice one can include other terms as well to speed up the convergence of the 
iterative routine. 

Now we want to prove the following statement: 

Theorem 1. If 5S[<j)(x)] = 0, then 4>(x) is a local minimum of the functional S[4>(x)}. 

Proof. By construction, the difference S[(f)] — S[<p] — oi\ | Ai |> 0, with equality possible only 
if Ai[(/>] = 0. Consider some family of functions 4>\(x), < A < 2 such that <fi\(x)\ x=1 = 4>(x). 
Then 

±S[fa(x)) |a=i=« 1 ^|A 1 [^(x)]| (27) 

because j^S[(j)\(x)] = at A = 1. 
Consider the following two cases. 

1) Suppose ^A 1 [0W(a:)] + 0. Since \k x \<^(x)\\ > 0, \k x \${x)\\ = inf ^[(j)^ {x)}\ = 0. 
Therefore, if ^A^^ (x))\ ^ 0, then (p(x) is a local minimum of S[cj)^ x \x)] with respect to 
A. 

2) If, on the other hand, ^|A 1 [0( A )(x)]| A =i = 0, then 



= _ Al [0W(x)] A=1 = -(T[^\x)]+2V[^ x \x)]) | A=1 = 

^(S[<P W (x)} + V[^\x)\) | A=1 = ^V[^\x)] x=1 (28) 

or, in other words, Ai can have a zero derivative only in the direction of the constant V. 
However, it was shown in Ref. J7| that the bounce <f>(x) is the global minimum of the action 
S[<f>] over the surface of constant V. 

In other words, the variation (-^(fi( x \x))d\ is tangential to the surface V = const, 
while the only 0] eigenvector of the second derivative operator (5 2 S/5<j)i5<f)j) = 5ijA + 
(d 2 U(x) I d(f>id(f)j) which corresponds to a negative eigenvalue was shown in Ref. to be 
orthogonal to the surface V = V[(f)] = const in the functional space. Therefore, is a 
minimum of S in this case as well. 
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We have shown that for an arbitrary variation ^ x '(x) near the bounce <f>(x), either (i) 
5S = 0, 5|Ai| p > 0, or (ii) 5S = SAi = 0, but 5 2 S/5(p 2 is positive definite. In either case, <p 
is the minimum of S. This completes the proof. □ 

Therefore, the improved action (|26|), as well as a more general form fl25|) has the desired 
solution 4>(x) as its minimum. This justifies the method described above of searching for the 
bounce as a minimum of the improved action rather than a saddle point of the actual action. 

Our method can be easily generalized to the case of finite-temperature field theory. In 
the high-temperature limit, the transition probability ||, 



3/2 



-S 3 [4>,T]/hT 



det'[-dl + U"(4>,T)] 



-1/2 



x(l + 0(h)), (29) 



det[-02 + tf"(O,T)] 

is determined by the three-dimensional action 5*3(0, T) of the 0(3)-symmetric solution <f>{x) 
of the equation 5S 3 = 0. As before, <j>(x) is a saddle point of S 3 [(f>}. The improved action 
53 [(f)], for which (j){x) is a minimum, can be constructed by analogy with equation (f25p: 



where 



SM := S z [<j>] + EKnIAfV + a 2 , n |Af ] r + a 3 ,^|AP J|p " 



Pn, a n > 



(30) 



AS T) = T + 3V 



(31) 



AP 



A 



{T)ij 



n 








/ <i 3 x det 






4>i 4>j 



U(4>i, ...,<f> n ,T) \<f > ( x )=<j > ( x ) 



(32) 



(33) 



The application of our method to the issue of stability of the charge and color conserving 
vacuum in the MSSM will be presented in II. 

The author would like to thank P. Langacker, A. Linde, G. Segre and P. Steinhardt 
for stimulating discussions. This work was supported by the U. S. Department of Energy 
Contract No. DE-AC02-76-ERO-3071. 
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